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We report the first observation of a B meson decay that is not accessible by a direct spectator 
process. The channel B° -> D+K~ is found in a sample of 85 x 10 6 BB events, collected with 
the Belle detector at KEKB, with a branching fraction B{B° -> D+K~) = (4.6±^ ± 1.3) x 10" 5 . 
We also obtain evidence for the B° — > D + ir~ decay with branching fraction B(B a — > D + ir~) — 
{2.4+1° ± °- 7 ) x 10_5 • This value 

may be used to extract a model-dependent value of \V u b\- 



PACS numbers: 13.25.Hw, 14.40.Nd 



Although the B mesons decay primarily through the 
spectator processes, other processes such as IF-exchange 
and final state interactions (FSI) may contribute appre- 
ciably, especially 0] for rare modes that are now emerg- 
ing at the B factory experiments. As the focus of CP 
asymmetry studies shifts to rarer modes, it is important 
to quantify the effects due to non-spectator and FSI pro- 
cesses, as they are often necessary for generating mea- 
surable CP asymmetries. Since the non-spectator effects 
carry significant theoretical uncertainties, a quantitative 
understanding must be developed before we can extract 
CKM matrix parameters. 

Important insight in this regard can be gained by ex- 
perimental measurement of channels where non-spectator 
processes dominate. We report here the first observa- 
tion of a mode that is not directly accessible through the 
spectator process. The mode B° — ► D + K~ may occur 
via W-exchange or final state rescattering, and predic- 
tions for its branching fraction vary over a wide range, 
3 x 1(T 6 - 1CT 4 |, g, |, §. The search for B° -> D+R- 
also encompasses B" — ► D + n~~ , a mode that is expected 
to be dominated by a (spectator) b — > u transition. As it 
lacks a penguin contribution, it can in principle provide 



a way to determine | V u b | @ ■ 

The results reported here are based on a 78.7 fb _1 
data sample, collected with the Belle detector || at the 
KEKB asymmetric energy e + e~ collider at the center- 
of-mass (CM) energy of the T(45) resonance and con- 
taining 85.0 x 10 6 produced BB pairs. A 7.5 fb _1 data 
sample taken at a CM energy that is 60 MeV below the 
T(45) resonance is used for systematic studies of the 
e + e~ — > qq background. 

The Belle detector has been described elsewhere || . 
Charged tracks are selected with requirements based on 
the average hit residual and impact parameter relative 
to the interaction point (IP). We also require that the 
transverse momentum of the tracks be greater than 0.1 
GeV/c in order to reduce the low momentum combi- 
natorial background. For charged particle identification 
(PID) the combined information from specific ionization 
in the central drift chamber (dE/dx), time-of-flight scin- 
tillation counters (TOF) and aerogel Cerenkov counters 
(ACC) is used. At large momenta (> 2.5 GeV/c) only 
the ACC and dE/dx are used. Charged kaons are se- 
lected with PID criteria that have an efficiency of 88%, 
a pion misidentification probability of 8%, and negligible 
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contamination from protons. The criteria for charged pi- 
ons have an efficiency of 89% and a kaon misidentification 
probability of 9%. All tracks that are positively identified 
as electrons are rejected. 

Neutral kaons are reconstructed via the decay K$ — ► 
7r + 7r~. The two-pion invariant mass is required to be 
within 6 MeV/c 2 (~ 2.5a) of the nominal K° mass, and 
the displacement of the tt + tt~ vertex from the IP in the 
transverse r-<f> plane is required to be between 0.1 cm and 
20 cm. The directions in the r-<j> projection of the K 2 
candidate's flight path and momentum are required to 
agree within 0.2 radians. 

We reconstruct D + mesons in the channels D + — ► 
07T+, K*°K + , and K$K + (inclusion of charge conju- 
gate states is implicit throughout this report). <f> (K*°) 
mesons are formed from the K + K~ (K + tt~ ) pairs 
with invariant mass within 10 MeV/c 2 (50 MeV/c 2 ) of 
the nominal <fi (K*°) mass. We select D + mesons in 
a wide (±0.5 GeV/c 2 ) window, for subsequent stud- 
ies; the M(D S ) signal region is defined to be within 
12 MeV/c 2 (~ 2.5cr) of the nominal D+ mass. D+ 
candidates are combined with a charged kaon or pion 
to form a B meson. Candidate events are identified 
by their CM energy difference, AE = E t ) — E b , and 
the beam constrained mass, Mb c = \J E\ — QZjPi) 2 , 
where Eb — y/s/2 is the beam energy and pi and Ei 
are the momenta and energies of the decay products of 
the B meson in the CM frame. We select events with 
Mbc > 5.2 GeV/c 2 and \AE\ < 0.2 GeV and define the B 
signal region to be 5.272 GeV/c 2 < M hc < 5.288 GeV/c 2 
and | AS | < 0.03 GeV. The M bc sideband is defined as 
5.20 GeV/c 2 < M hc < 5.26 GeV/c 2 . We use a Monte 
Carlo (MC) simulation to determine the efficiency ||. 

To suppress the large combinatorial background that 
is dominated by the two-jet-like e + e~ — * qq continuum 
process, variables that characterize the event topology 
are used. We require | cos6> t hr| < 0.80, where 6* thr is the 
angle between the thrust axis of the B candidate and 
that of the rest of the event. This requirement elimi- 
nates 77% of the continuum background and retains 78% 
of the signal events. We also define a Fisher discriminant, 
J 7 , that includes the production angle of the B candidate, 
the angle of the B candidate thrust axis with respect to 
the beam axis, and nine parameters that characterize the 
momentum flow in the event relative to the B candidate 
thrust axis in the CM frame ||. We impose a require- 
ment on T that rejects 50% of the remaining continuum 
background and retains 92% of the signal. 

We also consider possible backgrounds from qq events 
containing real D + mesons. These events peak in the 
M(D S ) spectra but not in the AE and M bc distributions. 
We study this background using the M bc sideband and 
find it to be fewer than 0.1 and 0.5 events for the B° — > 
D+K~ and B° — > D + tt~ modes, respectively. 

Other B decays, such as B° — > D + ir~, D + — ► 
K~it + it + , with one pion misidentified as a kaon, require 
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FIG. 1: The AE versus M bc scatter plot for the B° -> D+K~ 
(top) and B° — » D^ty~ (bottom) candidates in the M(D S ) 
signal region. The points represent the experimental data 
and the boxes show the B meson signal region. 
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FIG. 2: The M(D S ) spectra for B° -» D+K' (top) and 
B° — » D + w~ (bottom) in the B signal region. The points with 
errors represent experimental data and the curves display the 
results of the simultaneous fit described in the text. 



particular attention because they have large branching 
fractions and can peak in the M bc signal region. The re- 
constructed invariant mass spectra for these events over- 
lap with the signal D + mass region, while their AE distri- 
bution is shifted by about 50 MeV/c 2 . To suppress this 
background, we exclude event candidates that are con- 
sistent with the D + — > K~tt + it + mass hypothesis within 
15 MeV/c 2 (~ 3c) when the two same-sign particles are 
considered to be pions, independently of their PID infor- 
mation. For the D + — > KgK + mode there is a similar 
background from B° — > D + ir~, D + — ► K < gH + . In this 
case we exclude candidates consistent within 20 MeV/c 2 
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M^ K - for the 



K%K+. 



(~ 3cr) with the Z? + — > Ag7r + hypothesis. 

Possible backgrounds from i? decays via & — > c transi- 
tions (_B — > D S DX) are also considered. The from 
these decays have a lower momentum and are kinemati- 
cally separated from the signal. We analyzed a MC sam- 
ple of generic BB events corresponding to about twice 
the data sample and found no peaking backgrounds. 

Another potential BB background is charmless B° — ► 
K~K+K-tt+{K^K+K~). Such events peak in the AE 
and Mbc spectra, but not in the M(D S ) distributions. 
They tend to be dominated by quasi-two-body decay 
channels such as 4>K^° @ y|- 

To reduce this back- 
ground, we reject events with low (< 2 GeV/c 2 ) two 
particle invariant masses: M K - n + and 
Df — ► channel, M K + K - and M^»o K - for I? 
K*°K+, and M K+K - and M^o*- for L>+ 
The remaining background from these sources, if any, is 
excluded by fitting the M(D S ) distribution. 

The scatter plots in AE and M bc for the B° -> L>+ A~ 
and 5° — » D+ir~ candidates in the M(D S ) signal region 
are shown in Fig. [l]; a significant enhancement in the 
B signal region is observed. Figure || shows the M(D S ) 
spectra for selected B° -> D+K~ and B° -> L>+ 7r can- 
didates in the B signal region. In addition to clear signals 
at the Df mass in Fig. [2[ we observe peaks at the D + 
mass, corresponding to B° — > D + ir~ and B° — > D + K~ , 
D+ -> (fm + ,K*°K + ,K%K + . 

Our studies have shown that the backgrounds may 
peak in the signal region of M(D S ) or of AE (and M^ c ) 
but not in both simultaneously. To extract our signal, 
we therefore perform a binned maximum likelihood fit to 
the two-dimensional distribution of data in M(D S ) and 
AE, separating the backgrounds from the signal com- 
ponent, which peaks in both. For each of the three 
.D+ decay channels the AE range, —0.1 GeV< AE < 
0.2 GeV, is divided into 30 bins and the M(D S ) range, 
1.5 < M(D a ) < 2.5 GeV/c 2 , into 200 bins. All bins in all 
Df submodes are fitted simultaneously to a sum of sig- 
nal and background shapes. The signal is described 
by a two-dimensional Gaussian, with widths in both di- 
mensions obtained and fixed using reconstructed signals 
in the data from B° -> D+w-(D+ -> A""7r + 7r+, K s tt+). 
The signal amplitude is constrained to correspond to the 
same branching fraction B(B° — > D+h~) for all three 
Df submodes. The fit also includes an additional two- 
dimensional Gaussian for B° — > D + h~ decays. 

The background includes three components: combi- 
natorial (flat in M(D S ) and AE), qq events that peak 
in M(D S ) and are flat in AE, and B decays that peak 
in AE and are flat in M(D S ). The levels of the three 
components are allowed to vary independently in the 
three reconstructed Df modes. The fit results are given 
in Table |. The statistical significance quoted in Ta- 
ble | is defined as y/—2 \u(Cq / C max ) , where C ma x and 
Co denote the maximum likelihood with the fitted signal 
yield and with the signal yield fixed to zero, respectively. 




FIG. 3: The AE spectra for B° -> D+K~ (top) and B° -> 
Dfir~ (bottom) in the B signal region. The points with errors 
are experimental data and the curves are the results of the 
simultaneous fit described in the text. 



The results of one-dimensional fits to the M(D S ) and 
AE distributions are also shown in Table [I] for compar- 
ison. Figures || and |^ show the M(D S ) and AE pro- 
jections for events from the signal region and the fitted 
signal plus background combined shape by solid lines and 
background shape including the peaking background by 
dashed lines. The peaking background is found to be 
1.0 ± 0.5 and 1.6 ± 1.0 events for B° -> D+K~ and 
B° — * Z)+7r~, respectively. 

B° — > D* + h~ final states, where the low energy photon 
from the D* — > D s ^ decay is missed, can populate the 
B° — > D+h~ signal region. These would produce a long 
tail on the negative side of the AE distribution. In theo- 
retical models based on factorization, the B° — > £>* + 7r~ 



and B° — > Dfir decay widths are predicted to be ap- 
proximately equal; there are, however, no corresponding 
predictions for B v — * D s K~ decays. To study the 
sensitivity of the measured branching fraction to a pos- 
sible B° — * D* + h~ contribution, we perform a fit with 
an additional B — ► D* + h~ component included, where 
the signal shape is fixed from the MC and the branch- 
ing fraction is left as a free parameter. The resulting 
2% difference in the B° — > Dfir~ event yield (compared 
to the results presented in Table |) is added to the sys- 
tematic uncertainty; the change in the B° — ► D+K~ 
yield is less than 1%. We also check for crossfeed be- 
tween B° — ► DfK~ and B° — > Dfir~ due to kaon/pion 
misidentification. To study this we include the crossfeed 
contributions in the simultaneous fit, with shapes fixed 
from the MC and misidentification probabilities obtained 
from data; the uncertainty due to this effect is found to 
be negligible (< 1%). 

As a check, we apply the same procedure to B° — > 
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TABLE I: Results on the signal yields and branching fractions. The efficiencies do not include intermediate branching fractions. 



Mode 








M(D S ) - AE yield 


M{D S ) yield 


AE yield 


Efficiency, % 


B (10~ 5 ) 


Stat, significance 


go _^ 

B° -> 
B° -* 


u s I\ 

D+K- 
D+K- 


~,D+- 
~,D+- 
~,D+- 


— ► <jm~T~ 
-+ K*°K+ 
^K»K+ 


q+3.3 

2 7 

1 6+ 20 


a n+3-4 

».y_ 2 7 
5.1+1-5 

i-o±i:g 


q n+3.3 

2 7 

5.9^;° 


ii « _i_ n a 
6.8 ±0.3 
7.0 ±0.3 


- 1+2.0 , , 4 

4.8±^;| ± 1.3 
2.2+. 2 ;* ±0.6 


5 lcr 
3.8cr 

1.6(7 


B° -> 


D+K' 


, simultaneous fit 


16 4+ 4 ' 6 


15.0±t 6 8 


17.5±$» 




4.6+.1; 2 ± 13 


6.4cr 


B" -> 


D + 77- 


,D+- 




4 -'-2.0 


4.8±?;« 


4 n+ 2 - 6 

4 - u -2.0 


12.9 ±0.4 


2.4 + i;q ±0.7 


3.2cr 


B° -* 
B a -> 


D + 7T- 
D + TT- 


,D+- 
, D+ - 


4 K*°K+ 


3 4+ 32 

°- -2.4 

i-6±?:g 


2 Q+2- 8 
z - a -2.0 
o o+2.2 


44 +3.3 

u ' M -0.9 


7.5 ±0.3 
7.2 ±0.3 


2.4t 2 ;3 ±0.7 
2.2j;|;2 ±0.6 


1.6cr 

0.9CT 


B" -> 


D + 7T- 


, simultaneous fit 


io.it 4 ; 4 


io.3±*:i 


"-•^ 




2.4tg;° ±0.7 


3.6cr 



D+tt- and B° _-> D+X", L>+ -> 0tt+, +, 
and obtain -> D+tt-) = (2.8 ± 0.2) x 10~ 3 and 

#(B° -> D+K-) = (3.0_± 0.7) x 10" 4 , which agree well 
with the PDG values B(B° -► D+tt") = (3.0±0.4) x 10~ 3 
and 6(5° -► D+K-) = (2.0 ± 0.6) x 10~ 4 ||. 

The following sources of systematic error are found 
to be the most significant: tracking efficiency (2% per 
track), charged hadron identification efficiency (2% per 
particle), Kg reconstruction efficiency (6%), signal-shape 
parameterization (5%) and MC statistics (3%). The 
tracking efficiency error is estimated using r\ decays to 77 
and tt+tt^tt . The K/ir identification uncertainty is de- 
termined from D*+ — ► D°n+, D° — > K~it+ decays. We 
assume equal production of B+B~ and B°B° pairs but 
do not include an additional error from this assumption. 
The uncertainty in the D+ meson branching fractions, 
which is dominated by the 25% error on B{D+ — ► 4>tt+) 1 
is also taken into account. The overall systematic uncer- 
tainty is 28%. 

In summary, we report the first observation of B — ► 
D+K~ with a 6.4er statistical significance. We find 
B{B° -► D+K-) = (4.6 + ^ ± 1.3) x 10" 5 , which is 
consistent with a calculation of the IF-exchange rate 
in the "PQCD factorization" approach ||, but much 
higher than an earlier result Q. On the other hand, 
it should be noted that the recent observation of highcr- 
than-predicted rates for B a -> D°h° (h° = tt°, r), u) jlf, 
(i~§ Q and 5° -> L>°p° p) suggest that FSI may 
contribute appreciably to B° — > D+K~ fl], §. We 
also obtain evidence for £?° — > D+ir~ with £>(£>° — > 
= ( 2 - 4 ta8 + °- 7 ) x 10_5 ( 3 - 6cr statistical signif- 
icance). Our results are consistent with recent evidence 
from BaBar @. 

Since the dominant systematic uncertainty on both 
measurements is due to the branching fraction of D+ — > 
071+ , S^, we also report £?(£?° — > D+K~) x = 
(16.4±|;| ± 2.1) x 10- 7 and B{B° -> .D+tt-) x £ 07r = 
(8.6t^ ± 1.1) x 10~ 7 . Using B(B° -► D+tt~) /B(B° -> 
D+D-) - (0.424± 0.041) x |Kb/U cb | 2 @, and S(B° -> 
D+D - ) x S 07r = (3.0 ± 1.1) x 10" 4 calculated from a 
CLEO result JlSj, we can extract a model-dependent 
value iKh/Vctl = (8.2 + |^±3.4) x 10 -2 , where no error on 



the factorization assumption or other sources of model- 
dependence is included. This value is in agreement with 
ones obtained from inclusive semileptonic decays fl2|. 
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